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to turbulence ( weak ) immediately above tlie convection thresliold. 

Since the hydrodynamic equations involve coupled 
non-linear partial differential equations which are difficult to 
solve, a Lorenz-like truncation of the actual hydrodynamic 
equations is carried out which results in a system of coupled 
non-linear ordinary differential equations. This system provides 
an adequate description of convection near the threshold and 
facilitates a detailed study of the various convection states near 
the onset. 

In the present work, velocity, vorticity and the 
temperature fluctuations have been represented by Fourier modes 
upto second order. The critical Taylor number in the limit of very 

4 

large Prandtl number, has been determined to be 23-45 n by rolls 
at an angle of 58°with respect to the original set. This is in 
complete agreement with the results of Kuppers and Lortz. 

If the layer of fluid is assumed to be of finite Prandtl 
number, the non-linear advective term in the hydrodynamic equation 
cannot be neglected. Our truncated system makes it possible to 
arrive at closed form expressions even at low Prandtl numbers. The 
threshold Taylor number is found to be considerably reduced. This 
is in agreement with the numerical works of KOppers and the 
experimental works of Niemella and Donnelly. 

He then extend our study to consider the Kuppers Lortz 
instability for an electrically conducting fluid which is placed 
in an external magnetic field. T'he ooint effect of rotation and 
magnetic field on the onset of convection is particularly 



Kuppers-Lortz instability in this case shows the growth of a set 
of rolls at right angles to the original set at very low Taylor 
nunber indicating presunably the fornation of square rolls rather 
than cylindrical ones. 



CHAPTER I 


INTRODUCTION 

Instability in rotating fluids has drawn a lot of attention 
ever since the experimental evidence^ of the onset of time 
dependent motion immediately above the threshold for convection 

2 3 

was reported. Ktippers and Lortz and later Busse and Clever have 
analyzed that this direct transition to the chaotic state is 
related to the onset of an instability which is characteristic of 
rotating fluids. In this dissertation we investigate the Kuppers 
Lortz instability in a rotating Rayleigh Benard system using 
truncated models. 

Section I of this chapter reviews the concepts related to 
Kuppers Lortz instability and turbulence while Section II is a 
brief description of the work presented in the thesis. 

SECTION I 

REVIEW 


1 . 1 Turbulence 

A class of phenomena which is extremely sensitive to the 
variations in the initial conditions and is characterized by 

4 > 

randomness and irregularities , is termed as turbulence . For 

15 

fluids in particular .turbulence is related to a time dependent 
complex flow in contrast to a laminar flow in 


which the fluid 
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travels along a well defined path with a velocity which is 
spatially and temporally uniform. Beyond a critical Reynolds 
number the laminar flow becomes turbulent. After turbulence sets 
in, two flows which are identical with slightly different initial 
conditions may ultimately evolve into two absolutely different 
flows. Motion which is initially periodic may become 
quasi-periodic or aperiodic and finally acquire a highly random 
character . 

An extremely simplified explanation of turbulence is the 
absence of a unique solution of the equations of motion for 
asymptotically large times. When a system has stable solutions, 
turbulence does not set in -as the system settles into one of these 
stable states eventually. However, if all the solutions are 
unstable, the system is forced to migrate in phase space spanned 
by these solutions without settling into any of them. Thus a 
conceptual model of turbulence*^ is a manifold of stationary 
solutions all of which are unstable to some other solution in the 
manifold such that the system constantly moves from the 
neighbourhood of one solution to another but never manages to 
settle into any of them. Therefore the phenomenon appears to be 
completely dependent on the qualitative features of the 
differential equations governing the system. The basic observation 
is that simple systems"*^ of three or more coupled non-linear first 
order equations often admit chaotic solutions called strange 
attractors. It is this attractor or a few attractors embedded in 
the state space of the fluid which account for the turbulent 
behaviour. Motion on these attractors are highly sensitive to 
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initial conditions and possibly this sensitivity accounts for the 
stochastic time dependence of the fluid and the unpredictability 
of the time evolution of the system. 

Though turbulence is a very common phenomena in nature we 
shall be concentrating only on simple physical systems under 
controlled conditions for the study of onset of turbulence. The 
dynamical behaviour of some physical systems is dependent strongly 
on a tunable parameter called the control parameter (r). The 
system remains in a steady state for low values of r. At r = r , 
the critical value, this steady state may be replaced by another 
steady state or an oscillatory state. Increasing r further, may 
result in the onset of time dependent motion finally leading the 
system to a chaotic state. A simple laboratory arrangement for the 
study of onset of turbulence is the Rayleigh Benard system which 
is described below. 

1 . 2 Rayleigh Benard system 

A fluid of positive expansion coefficient, confined to a 
rectangular box of depth 'd'is considered. For theoretical 
simplifications, the lateral dimensions of the box are assumed to 
extend to infinity. The fluid is heated uniformly from below. The 
physical conditions of the fluid are described by two 
dimensionless parameters, the Rayleigh number R 

a g d^(T^-Tp 


R 
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and the Prandtl number o' = 

« 

where ot = coefficient of volume expansion 
g = acceleration due to gravity 
« = coefficient of heat diffusivity 
V = coefficient of kinematic viscosity 
= temperature of lower plate 

= temperature of upper plate 

Since the fluid is being heated from below, an adverse 

temperature gradient exists along the layer which results in a top 

heavy arrangement . This is a potentially unstable configuration. 

The natural tendency of the fluid to redistribute itself and 

attain stability is opposed by its viscosity and conductivity. 

This stationary state of the fluid however gets destabilized at a 

critical temperature gradient, characterized by the Rayleigh 

number R . The energy released by the buoyancy force overcomes the 
c 

dissipation due to the viscous forces and the fluid sets in motion 

with the hot fluid rising and cold heavy fluid flowing down. Thus 

for the Rayleigh Benard system at r = we have the onset of the 

first instability where the stationary conduction state 

destabilizes against steady state convection. Experimentally it 

has been observed^ that the steady state convection state remains 

stable till about r = 12 r . Above this value of r time dependent 

c 

motion sets in finally leading the system to a chaotic state. To 

determine this r we need to first concentrate on the equations 

c 

guiding the motion of the fluid. 
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1.3. The hydrodynamic equation in the Boussinesq approximation 

As the convecting fluid moves across the layer, its physical 
properties like the viscosity and density are likekly to change 
with the variation in temperature. In the Boussinesq 
approximation^, variations of all fluid properties except the 
density are ignored completely. Variations in density are also 
ignored except in so far as they give rise to the buoyancy force. 
The conservation of mass, momentum and heat are expressed 

respectively as 

-> -*■ 

the continuity equation 7-V = 0 (1.1) 

the Navier Stokes equation 

->• 

^v ->->-> VP-^ <5p 

+(v V>v = g ( 1+ — — ) +v V ^ V (1.2) 

d t P 

o 

where <^p is the variation of density from the equilibrium point. 
The equation of heat diffusion in the temperature variable T is 

^ T ^ 

+ (v.V )T = « V ^ T (1.3) 

9 t 

These equations can be reduced to a dimensionless form by 
scaling 1) all distances by d 

2) velocity by «/d 

3) time by d^/« 

4) temperature by AT (AT=T 2 -T^) 

We describe the Rayleigh Benard convection in terms of the z 
component of velocity (w) and O, the temperature fluctuation from 
the equilibrium state. If we consider fluids of very high Prandtl 
number, the non-linear term in the Navier Stokes equation can be 
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neglected and we arrive at the dimensionless form of the 
hydrodynamic equations (Appendix A) 


V 4 , 


- 7 


d w 


e t 


-) = 


-R 


(l-4a) 


V 


2 


e 


a 


d e 

9 t 


-w + (vV ) 0 


(1.4b) 


V. r, 2 9 ^ , 0 ^ 

where v = + 

9 rT 9 Y 

The solutions of equations (1.4a) and (1.4b) have to satisfy 
certain boundary conditions depending on the type of bounding 
surfaces within which the fluid is confined. Irrespective of the 
type of boundary we have 

w = 0 at z = 0 and d (1.5a) 

If the boundaries are rigid we have the no slip condition 
on the velocity field. This is equivalent to having the x 
component (u) and the y component (v) of the velocity field to be 
zero on the boundaries. 

u=0,v=0 (1.5b) 

Since this condition has to be satisfied for all x and y on the 


surface we have — 5 = 0 

9 2 

For stress free boundaries the vanishing of 

stress would imply — ^ 1 ^ — = 0 

9 -z 9 z 


and using the continuity equation we have 



(1.5c) 
the tangential 
(1.5d) 

= 0 (1.5e) 


If the bounding surface is a perfect conductor then the 
temperature fluctuation O has to be zero on the boundaries. 
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0 = 0 

However, if the boundaries are considered to be 


(1.5f) 

thermal 


insulators, we have — 3 = 0 

a z 

The equilibrium state of the system is described by 

-> 

V = 0 


(1.5g) 


(1.6a) 


V^T=0 (1.6b) 

T - 

T = ^ (l-6c) 

The stability of the equilibrium state (conduction state) has 
to be tested against infinitesimally small perturbations. For this 
purpose we employ the technique of linear stability analysis. 


1.4 Linear Stability Analysis 


In linear stability analysis, the effect of small 

7 

fluctuations away from a solution to the equations is examined as 
a function of a parameter (R, the Rayleigh number in our case). 
The hydrodynamic equations are linearized in the amplitude of the 
fluctuations and these linearized equations are solved to 
determine whether the perturbations decay or grow in time. 
Consequently, the boundaries of a stable region in the parameter 
space can be defined. 

Convective state 

Performing linear stability analysis on the conduction state 
determines the critical Rayleigh number for the destabilization of 
this stationary state. This takes place via a marginal state which 



8 


could be either tine independent or oscillatory. For a Rayleigh 
Benard system with no external constraint, the onset of 

B 

instability is always via the stationary state 

For the laterally infinite system that is considered, the 

convection pattern can be assumed to be periodic in the x-y plane 

with a wave number a . The critical wave number a is obtained by 

c 

finding the minimum of R(a). The critical Rayleigh number is 

7 

obtained by setting a = in the expression determined for the 

critical Rayleigh number. 

Linear stability analysis determines the magnitude of the 

wave number. It provides no information about the horizontal 

structure of the convective motion^. This critical wave number 

that is determined can be resolved in an infinite number of x and 

y components and corresponding to each there is a steady state 

solution. In order to determine the form and stability of the new 

steady state solutions which evolve from the instability, the 

non-linear equations have to be considered. They are expanded in a 

1 1 

sequence of inhomogeneous equations . From the solvability 
criteria described in Appendix B, solutions of higher orders are 
arrived at. 

Employing the perturbation technique, Schldter, Lortz and 
12 

Busse have concluded that the inclusion of the non-linear terms 

in the analysis lifts the degeneracy in a completely. They have 

c 

? 

also concluded that the only physically realized form of 
convection near the threshold is in the form of steady state two 
dimensional rolls (fig.l). All other forms of the solution of the 
non-linear equations decay exponentially with time. As mentioned 
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earlier, the system remains in the steady state over a large range 
of Rayleigh number and does not exhibit any characteristics of 
turbulence . 

However, if the Rayleigh Benard system is subjected to an 
external constraint, as for example, rotation about a vertical 
axis; the stability picture undergoes a complete change. For R 
slightly greater than R^, irregular time dependent motion has been 

13 

observed . This system thus appears to be a simple model 

3 

exhibiting the characteristics of a weakly turbulent system . 

In addition to the Rayleigh number and the Prandtl number, 
the fluid is also described by another dimensionless parameter, 

4 2 2 

the Taylor number T = 4 d Ci / x> where O is the angular velocity 
of the fluid and the other symbols have been defined in the 
previous sections. 

The Navier Stokes equation has to be augmented by the terms 
for the Coriolis and centrifugal forces while the equation for 
heat diffusion remains unaltered. The equivalent form of Navier 
Stokes equation is 


d V 
d t 



1 

2 


-> -> 
|0 X r 


2 


) + V V 


2 


V +2 V X Q 


- ( 1 + 



(1.7) 


Using the same scaling factors as for the non-rotating case we 


14 


arrive at the equations (Appendix A) 


2 2 
V (V 


a t 


) w = 


a c 


a z 


R 7 ^ 0 
1 


(1.8a) 
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2 


d 

d t 


) C 


( 7 


2 


a 


d 

d t 


) 0 = 


where = T 


1.5 Conductive state 


-T 


d w 
d z 


- w + ( v.v ) e 


(1.8b) 

(1.8c) 


The equilibrium state is defined by 
V = 0 , T = T 

S 

In order to study the stability of this conduction state as 
in the non-rotating case, we perturb this equilibrium state 
infinitesimally. The growth of these perturbations are determined 
by solving the linearized hydrodynamic equations in these 
fluctuations. An exponential time dependence e^^ is assumed for 
the perturbations. The conduction state would destabilize against 
these perturbations if Re (p) > 0 . Unlike the non rotating case, 
the marginal state can be either stationary or oscillatory 
depending on whether Im (p) = 0 or not for Re (p) = 0. For Im (p) 
not zero, the instability grows with an oscillating amplitude. 
However if the analysis is confined to fluids of large Prandtl 
numbers, the onset of instability is always stationary. 
Oscillatory instability is always preferred for fluids with 
Prandtl number less than 0.67 For fluids with high Prandtl 
number since the oscillatory state is ruled out in the non-linear 

IB 

region, we have a direct transition to a chaotic state . 

Linear stability analysis predicts that the conduction state 
is destabilized in favour of steady state convection at a 
critical Rayleigh number which is dependent on the Taylor number. 
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7 

Rotation thus has an inhibiting effect on the onset of convection 

implying that any finite value of T increases R beyond that 

calculated for the non rotating case. 

As in the non rotating case, linear stability analysis 

predicts just the magnitude of the critical wave number. Due to 

the horizontal isotropy of the system all directions are equally 

preferred. Thus an infinite number of solutions all with the same 

wave number but differing in their orientation is possible. For 

determining the physically realized solution, inclusion of the 

non-linear terms in the equations is essential. Employing the 

1 1 

perturbation technique, SchlUter, Lortz and Busse and Kuppers 
2 

and Lortz have determined that two dimensional rolls represent 

the only form of stable solutions for the rotating case as well. 

However, unlike the non rotating case, these steady state two 

dimensional rolls become unstable if the Taylor number exceeds a 

critical value T . This happens for R just above the threshold. 

c 

An investigation of this instability of the convection rolls for 
slightly super-critical Rayleigh number is the key to studying the 
Kuppers Lortz instability which is discussed in the next section. 

1.6 Kuppers Lortz Instability 

Rayleigh Benard convection in a rotating horizontal layer 
develops an interesting complication just above the threshold. 
Kuppers and Lortz^ have discovered that a set of cylindrical rolls 
formed at the onset of convection destabilizes in favour of 
another set developing at an angle to the original roll provided 
the Taylor number is above the critical value T^. 


These 
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perturbing rolls do not grow indefinitely. They in turn sucuumb to 
a similar type of instability. Thus we have a situation in which 
the solutions switch over from one roll system to another^. This 
transition from one set to another is nearly periodic with a 
transition time that fluctuates statistically about a mean value 
thus giving rise to weak turbulence. 
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SECTION II 


QYERYIEH 


In the chapters that follow we investigate the Kiippers Lortz 
instability in horizontal fluids using truncated models. While 
the linearized hydrodynamic equations are adequate for studying 
the onset of convection, the non-linear equations are essential 
for tracing the course of the system after the onset. The coupled, 
non-linear partial differential equations are difficult to handle. 

Q 

With the discovery by Lorenz that the essential information is 
contained in just a few modes of the system, it has been possible 
to consider the non-linear effects in a simplified manner. Using a 

1 <4 

truncated representation one can reduce the partial differential 
equations to a system of coupled non-linear ordinary differential 
equations. The truncated system provides an adequate description 
of the onset of convection and can be numerically integrated 
easily (Appendix C) facilitating the study of non linear effects 
near the threshold. 

Chapter II deals with the Kuppers Lortz instability in fluids 
with infinite Prandtl number. We have first considered a minimal 
representation of the velocity, vorticity and temperature modes. 
This has the advantage of being handled analytically and hence the 
qualitative features of Kdppers Lortz instability can be explored 
in a simplified manner using such a model. However for more 
accurate results, we need to consider a more complete model by 
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considering higher order terms in the expansion of the velocity, 
vorticity and temperature fields. 

In chapter III we consider Kdppers Lortz instability in 
electrically conducting fluids in the presence of a magnetic 
field. Acting independently, both magnetic field and rotation 
have an inhibiting effect on the onset of convection . However, 
when the fluid is subjected to these two constraints 
simultaneously, the effect is to reduce the threshold Rayleigh 
number below the simple superposition of the two effects acting 
individually. Thus an interesting question is how would the 
magnetic field affect the Kdppers Lortz instability, tfe find that 
it always stabilizes the system against this instability. In fact 
the critical Taylor number is found to be a monotonically 
increasing function of the magnetic field. Presence of magnetic 
field also enables us to investigate Kuppers Lortz instability in 
low Prandtl number fluids. In the absence of magnetic field for 
Taylor number sufficiently high, the onset is oscillatory if the 
Prandtl number is less than 0.67. However if it is subjected to 
just the magnetic field the onset as stationary convection is 
ensured under all terrestrial conditions. If the magnetic field 
and rotation are present simultaneously there exists a range of 

magnetic field for which the onset is stationary' thus making it | 

[ 

possible to investigate Kdppers Lortz instability in low Prandtl 
number fluids particularly in liquid metals like Mercury (o' = 

0.025) I 

Finally in chapter IV we discuss the case of a rotating fluid 

i 

j 

layer between two thermally insulating boundaries. He determine [ 
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the critical Rayleigh number, the critical wave number and address 
the question of Kuppers Lortz instability. He find that for a 
range of rotation speeds, the onset of stationary convection is in 
the form of long wavelength rolls. An investigation of Kuppers 
Lortz instability shows that a set of rolls grow at right angles 
to the original set at very low Taylor numbers indicating 
presumably the preference of square rolls over the cylindrical 
rolls . 

As we increase the rotation speed, the rolls become unstable 
for a range of angles and finally for very large T the basic rolls 
would be unstable to rolls at all angles. 
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CHAPTER - II 


KOPPERS LORTZ IMSTABILITY IN FLUIDS WITH LARGE PRANDTL HUMBER 
2 . 1 The Lorenz model for Rayleigh Benard convection 


A layer of fluid of depth d and infinite lateral dimensions 
is heated from below such that an adverse temperature gradient is 
maintained across the layer. Though this represents an unstable 
configuration, the viscosity and conductivity of the fluid tend to 
stabilize the system and the fluid remains stationary. As the 
temperature gradient across the layer increases, the buoyancy 
force overcomes the dissipative effects of viscosity and motion 
sets in. This critical temperature gradient (defined by the 
Rayleigh number) also depends on the physical properties of the 


fluid. 


R = 


a g d'^ A T 


« V 


where ex is the coefficient of volume expansion 
g is the acceleration due to gravity 
AT is the temperature difference across the layer 
» is the coefficient of thermal diffusivity 
V is the coefficient of kinematic viscosity 
The motion of the convective fluid is guided by the Navier 


Stokes equation which in the Boussinesq approximation is 
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d V 


d t 


-> 

( V 


V ) V 


V P 


2 

V V 


( 1 + 


6 p 


) I 


( 2 . 1 ) 


where P represents the pressure, 6p is the variation of density 
from the steady conductive state. Considering T as the temperature 
variable, the equation for heat diffusion is 


d T 


d t 


(V 


-> 

V ) T 




2 

7 ^ T 


( 2 . 2 ) 


The equation of state is 


p^ { 1 + a (T^- T) } 


(2.3) 


Where P^is the density of the fluid at the lower plate which is at 
temperature T^. 

We assume the fluid to be surrounded by free surfaces and 
perfectly conducting boundaries. Though the assumption of the 
boundaries being free surfaces is unphysical, nevertheless it 
helps in determining the qualitative features of the flow in a 
simplified manner. Rigid boundaries are more likely to be 
encountered in the experiments and hence are more realistic but 
they have the disadvantage of involving a lot of mathematical 
complications. In most of our analysis we shall confine ourselves 
to free boundaries . In Appendix D we perform our calculations for 
the rigid boundaries. 

The conduction state of the fluid is described by a linear 
temperature profile and v = 0. When the fluid starts convecting 
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this linear temperature profile gets distorted due to the 
motion of the fluid > If 0 is assumed to be this deviation, the 
altered temperature is given as 


T = T 


T - T 
2 1 


9 


(2.4) 


With SP as the change in pressure, the equation of the 
convecting fluid becomes 


d V 


3 t 


( V 


V ) V 


■* 

- V ( 


<5 P 


) + g « 0 


+ V ^ V 


(2-5a) 


a 9 


d t 


2 

«v‘^e = (v.v)0 


w 


(2.5b) 


Performing the same scalings as in Chapter I, we arrive at 
the dimensionless form of equations (2.5a) and (2.5b). In 
Append ix-A we reduce these equations to a scalar form. Thus the 
system is described by w, the z component of velocity and Q, the 
distortion in the temperature field. The relevant equations are 


^ " ( 


a t 


) w 


-R V ^ 0 
1 


(2.6a) 


( V 


- a 


a 


a t 


) © 


-> 

( V 


-» 

V ) 0 


(2.6b) 


The solutions of these equations in general are represented in 
terms of the complete set of normal modes of the system. Lorenz 
truncation® refers to expressing these solutions in terms of only 
a few modes of the system. For the present problem these modes 
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are the Fourier modes. The simplest expansions for w and &, 
compatible with the boundary conditions (stress-free ) are 

w = a(t) cos a x sin rxz (2.7a) 


d = b(t) cos a X sin nz + c(t) sin 2nz (2.7b) 

Substitution of the expansions (2.7a) and (2.7b) in equations 
(2.6a) and (2.6b) results in a set of ordinary differential 
equations which are called the Lorenz equations. With X , Y and 
Z representing the scaled versions of a(t), b(t) and c(t) 

B 

respectively the dynamical system that is obtained is 

X = c?(-X + Y) (2.8a) 


Y= -XZ+rX-Y 


(2.8b) 


Z = X Y - b Z 


(2.8c) 


where b 


4 n 


/„2 ^ 2 . 
(fTf + a ) 


and r = 


R 


R 


We rescale time by a factor of o'. 

In these equations, X is proportional to the intensity of 
convective flow, Y represents the temperature difference between 
the ascending and descending currents. The variable Z is 
proportional to the distortion of the temperature profile from the 
linear distribution. 


2 . 2 Linear Stability Analysis 

The steady state solution of equations (2.8a)-(2.8c) is 
obtained by considering all time derivatives to be zero. In 
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general we consider the steady state solution to be represented by 
( determine the stability of this steady state we 
superimpose infinitesimal perturbations (x y z) on these 
solutions and study the time evolution of these perturbations by 
setting the hydrodynamic equations in them, neglecting all 
non-linear terms. If a time dependence of the form e^^ is 
assumed, the steady state would be destabilized by these 
disturbances if Re (p) is greater than zero. These perturbations 
are governed by the linearized equations which are 




X 

II 

Q 

1 

X + 

y ) 

(2.9a) 


y 

= ( 

r - ) X 

- 

y - z 

(2.9b) 



z 

= ^ 

4- 

X^ y - b z 

(2.9c) 

For a 

non 

trivial 

solution 

of 

these equations, the 

determinant 

of the 

matrix 

formed by 

the 

coefficient of 

x,y,z must 

necessarily 

be zero. Recalling that 

the 

perturbations 

have a time 

dependence i 

of the 

form e^ 

the solvability condition 

is given by 


the characteristic equation of the matrix 

(P + b){p^ + p(l + i^) + o'(l-r)} = 0 (2.10) 

with ( X Y Z ) = 0 for the conduction state, 

o o o 

For r > 0, the equation has three real roots. If r < 1 , all the 
roots are less than zero and the perturbations decay with time. 
For r > 1, one of these roots becomes positive and thus the 
conduction state gets destabilized at r = 1 ( R = For r > 1 

we have two additional steady state solutions given by 
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X = Y=: j:{b(r-l) 

Z = (r - 1) 

For both these solutions, the characteristic equation of the 
matrix is 
3 2 

P + p (b+l + o') + pb(r + ty) + 2b<7(r-l) = 0 

For r > 1 , this equation possesses one real negative root and two 
complex conjugate roots. 

If the real root is p^ and the complex conjugate roots are 
represented by + i we have 

(p-PQ)(p-ic*^Q)(p+iw^) = 0 (2.12a) 

or P^ - P^ P^ + P - P =0 (2.12b) 


Thus for the roots to be pure imaginary the constant term is the 

2 

product of the coefficients of p and p. Applying this condition 
to equation (2.11) we obtain the critical r for the onset of 
oscillatory instability as 


o-(b + 3 + cy) 
^os “ ( O' - b - 1 ) 


(2.13) 


Ifo’<(b+l) , the steady state convection rolls will always 
remain stable. For o' >(b+l)atr = r these rolls are 

, OS 

0 

destabilized and oscillatory instability sets in . 

2.3 Rayleigh Benard convection in a rotating fluid 


If the Rayleigh Benard system is rotated about a vertical 
axis , the region of stability gets considerably altered. Unlike 
the non rotating case, the marginal state can now be oscillatory 
provided the Prandtl number is less than 0.67. Another 
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characteristic feature of the rotating system is the appearance of 

time dependent motion immediately above the threshold. Kiippers and 
2 

Lortz have determined that among all the possible steady state 
solutions just above the threshold, only convection in the form of 
two dimensional rolls are stable. A similar observation was made 
by Schluter, Busse and Lortz^^ for the non rotating case. However 
for the rotating case, above a critical rotation speed 
characterized by the dimensionless parameter T, the Taylor number, 
these steady state two dimensional rolls are destabilized and no 
stable stationary solution is possible. In the limit of large 
Prandtl number, since oscillatory instabilities are precluded, the 
physically realized convection presents an enigmatic problem. 
According to Kuppers and Lortz, the primary cause of the 
restricted stability domain of the two dimensional convection 
arises from a non linear effect. The cylindrical rolls that are 

formed at the onset of convection get destabilized by another set 

of rolls which develop at an angle to the original set if the 
rotation speed is above the critical value. This instability which 
occurs right at the threshold is known as the Kuppers Lortz 
instability . For fluids with large Prandtl number the critical 

Taylor number determined was 2285 with the assumption of free 

boundaries. 

2.4 Kdppers Lortz instability 

An investigation of the time evolution of the KOppers Lortz 
instability gives an insight into the origin of the irregular 
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time dependent motion and the transition to a chaotic state above 
the threshold .We represent the steady state rolls by the 
expansions of the form^^. 


N 


w = f ( z , a ) 


i k 


. r 




t ) 


n =- N 


( 2.14 ) 


for the velocity field and similar expansions for the vorticity and 
temperature modes .Here a refers to the critical wave number. The 

C 

c^ ' s are assumed to be weakly time dependent . The horizontal 
vectors k ' s are arbitrary except that I k | = a 

^ ' n ' C , 

k = - k 

— n n . 

We start with the initial value for given by the steady 
state solution for the two dimensional rolls and much smaller 
values for all other coefficients ( fig. 2 )^. The coefficient 
corresponding to the vector k^ grows until it reaches an 
amplitude comparable to that of at which point , starts to 

decay and is eventually replaced by as the approximately steady 
state solution . At the same time , the disturbance amplitude 
corresponding to k^ starts growing exponentially . After a while 
the amplitude of C_ becomes comparable to that of C and finally 

3 ^ 

is replaced by C^. Since had been decaying before it started 
growing, it takes longer for this disturbance to reach the 
equilibrium value. This process is repeated cyclically but the 
period tends to grow without limit. This however is an unphysical 
feature of the solution. The properties of the system appear to 
depend on the time elapsed since the initial conditions were set. 
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This inhomogeneity in time, of the solution, can be resolved by 
realizing that the disturbances are present at all times and not 
just at the initial moment as assumed . The presence of noise 
level prevents the amplitudes from decaying to arbitrarily low 
levels. At the same time it introduces a random element into the 
time dependence of the system. Thus the weakly non-linear problem 
of convection in rotating system exhibits a typical feature of 
turbulent fluids in that its time evolution is stochastic in 
nature . 


2 . 5 Hydrodynamic equations for the rotating system 

Rayleigh Benard convection for rotating fluids is described 
in terms of three dimensionless variables w , the z component of 
velocity, C , the z component of vorticity and &,the deviation of 
the temperature from the conduction state. The motion of the 
convecting fluid is guided by the usual Navier Stokes equation 
augmented by the Coriolis and centrifugal forces . Thus we have 


d V 


a t 


+ ( V 


V ) V 


V 


IQ X 


r 


r, 2 

U V V 


2 V X O 


g ( 1 + 


6 p 


( 2.15a ) 


The equation for heat diffusion remains unaltered by the 
inclusion of the effects of rotation . Thus we have 
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d T 


a t 


( V 


V ) T 




2 

V ^ T 


( 2.15b ) 


As derived in appendix A , considering the curl on both sides of 
equation ( 2.15a ) and taking its z conponent results in the 
equation for C -We therefore have 


a C 


a t 


2 O 


a w 


a 2 


2 

V V C 


( 2.15c ) 


Using the results of Appendix A and the scaling factors of 
chapter I, with C scaled by « / d^, we arrive at the following 
dimensionless form of the equations. 


2 2 
V ^ ( V ^ 


a t 


) w = 


a c 
a z 


- R 0 (2-16a) 


(2. 16b) 

/ V cf — ^ = -w + (v.'7)Q (2.16c) 

a t 

We represent the steady state cylindrical rolls along the y 

axis in the form of 

w = a (t) cos ax sin nz 
C = f (t) cos ax cos 

Q = b (t) cos ax sin ^rz + c (t) sin 2nz 
Substitution of these expansions in equations (2 . 16a)-(2 . 16c) 
lead to a set of Lorenz equations similar to the non rotating 
case. Thus the dynamical system that we obtain is 


( ^ 


a t 


) C 
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<y i - X + Y 4 - (}t) 


Y = - XZ + rX-Y 


(2.17a) 

(2.17b) 


Z = X Y 


b Z 


G = -o'(G+tX) 


(2.17c) 

(2.17d) 


where b = 


4 


2 2 
+ a 


t2 . 


T 


(n^ + a^) ^ 


He scale time by . 

The steady state solution of this system of equations is given by 
X=Y=Z=G=0 .Here X, Y, Z and G are the scaled forms of a, 
b, c and f respectively. He perform linear stability analysis on 
this state to determine the critical Rayleigh number and the 
critical wave number for the onset of convection 
perturbations x,y,z and g with time dependence of the form e 
and realizing that the equation for z decouples when we linearize 
in these perturbations , we obtain the condition 

r = 1 + t^ 

for the onset of stationary convection . The critical wave number 

a is obtained by minimizing R(a) with respect to a. This leads 

c ; 


Assuming 
pt 


2 m^ + 3 m^ = 1 + 


(2.18) 


n 


to the condition 

h r 

where m = — j 

Substitution of a = a in the expression for R determines the 
critical Rayleigh number. 
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( )^ + T rr^ 

^ 1 ( 2 . 19 ) 

a 

c 

2 . 6 Kuppers Lortz instability with truncated model 

The cylindrical rolls that are formed at the onset are 
represented as 

= ( 1 + (2.20a) 

G = - t X (2,20b) 

o o 

= r - ( 1 + ) (2.20c) 

Investigation of the stability of these convection rolls 
forms the basis for the study of Kiippers Lortz instability. These 
rolls are stable against self perturbations, that is, by the rolls 
whose axis is along the y direction . We superimpose another set 
of rolls , the axis of which is in an arbitrary direction and 
study the time development of these rolls. The velocity, vorticity 
and temperature modes would then have y dependence as well. 
Accordingly we represent them by 

w = sin nz ( a(t) cos ax + a^(t) cos (k^x +k 2 y) } (2.21a) 

C = cos nz i f(t) cos ax + f^Ct) cos (k^x tk^y) } (2.21b) 

S — sin nz { b(t) cos ax + b^(t) cos (k^x + k^y) } 

+ sin 2nz C c(t) + c^(t) cos (k^x + a x + 

c (t) cos (k X + k y -ax)} (2.21c) 

z * ^ 
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The above expansions are consistent with the boundary 
conditions for stress free, perfectly conducting surfaces with 
2 

w = — - — ^ — = 0 , =0,0 = 0 (2.22) 

Introducing these expansions in equations (2.21a)-(2.21c) we 
obtain the corresponding Lorenz model. Scaling time by and 


4 + 2 a^(l ± cos &) 


,2^2 .0.5 
( n + a ) 


2 2 
( ^ + ) 


n 


we have 


X = <y(-X + Y + tG) 


(2.23a) 


Y = -XZ + rX-Y- 

Z_ 


+ c sin & 


c sin 0 


} 


I ( 1 - cos e ) x^ 
( 1 + cos 0 ) X^ 


} 


(2.23b) 


Z = X Y + X Y^ - b Z 
1 1 

G = - ©-(G + tX) 

= c/ ( - X^ + Y^ + t G^ ) 

Y^ = - X^ Z + r X^ - - 


(2.23c) 


(2.23d) 


(2.23e) 


4 


{ ( 1 - cos 0 ) X 
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- G c sin 9 } 



1 + cos 9 ) X + G c sin 9 } 


(2.23f ) 


( 1- cos 9 ) 


X + 


X. Y ) 


K -7 , c sin 9 - 

*>+ + 2 ' 


Y - G y^) 


(2.23*) 


L i _!-. sag g . ) < X Y, + X, Y ) - b_ ^ (G Y,-G^ Y) 


G^ = - O' ( X^ + G^ (2.23i) 

The two sets of rolls are given either by X,Y, Z and G non 
zero and the others zero or by X^ , Y^ , Z and G^ non zero and the 
others zero. We assume the basic pattern to be the rolls along the 
y axis which is given by the former conditions and the other set 
is assumed to be representing the perturbing set. Linearizing in 
these perturbations, we set up the hydrodynamic equations with the 
assumption that they have an exponential time dependence. The 
condition for the onset of Kiippers Lortz instability is given by 

P = 0. 

The rolls along an arbitrary direction are described by 


Y^ = ( 1 + t^ ) , G^ = - ^ \ (2.24a) 

After dropping all time derivatives, 
expressions for Z^ and Z^, as 


we obtain the 
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- Cl- cos ^ ) ( 1 + ) X (2.24b) 

b_ = ( 1 + cos 0 ) ( 1 + ) XX (2.24c) 

1 

The condition for KOppers Lortz instability is given by 
considering = 0 . This is equivalent to 

^^(1 - cos 0 + c t sin 0 ) + j^(l + cos 0 - c t sin 0 ) = 0 

(2.25) 

where we have used ( 2.20a)-(2.20c) and (2.24a). Substitution 

of (2.24b) and (2.24c) results in the determination of the domain 

of Kuppers Lortz instability . The critical Taylor number is 
obtained as^^ 

(1 + cos 0)^ b + (1 - cos 0)^ b_ 

c t sin 0 = (2.26) 

(1 + cos 0) b, - (1 - cos 0) b 

Substitution of c ,t , b^ and b_ results in the critical wave 

1 s> 

number as a as (V3 )n and corresponding to this we have T as 

c ^ 

80 rr*. Though this result is somewhat higher than that calculated 
by Ktippers and Lortz but the minimal model is useful in 
determining the qualitative features of Kdppers Lortz instability 
. The critical angle of the perturbing rolls predicted by this 
model is ei**. For more accurate results , we consider a higher 
order truncated model for the velocity and vorticity modes. 
Accordingly in w we have the additional terms of the form 
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sin 2nz {a^C t)cos(k^x + ax + k^y) + a^(t)cos(k^x - ax + k^y)} 
Similarly the additional terms in the expression for C are 
cos 2a2{f (t)cos(k X + ax + k y) + f (t)cos(k x - ax + k v)} 

^ ^ ^31 2 

The expansion for O remains unaltered and is given by equatioi 
(2.21c) . Inclusion of higher order terms in the expansion of w 
and C gives rise to a more complicated Lorenz model. Our 
dynamical system now contains a set of 13 ordinary differential 
equations. The additional equations are 


2 t G. 


1 “b 


2 Z^(l + cos ©) 


\ } 


(2.27a) 


2 t G. 


I ^ i. 

1 "b 


2 (1 


cos 0) 


(2.27b) 


G^ = - cy ( 2 t ^2 ^ (2.27c) 

G = - <y(2tX^ + bG, ) (2.27d) 

3 3-3 

MT » * * 

The equations for Y, Z^, Z^ are altered due to the 

additional terms in the velocity and the vorticity modes. 
However, the equations for X , G , X^, G^ and Z remain unaltered 
and are given respectively by (2.23a), (2.23d), (2.23e), (2.23i) 
and (2.23c). The altered equations are 
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Y = -XZ +rX-Y- 


4 { 


( 1 - cos © ) X 

1 . 


+ c sin © G, 


} 


{ 


4 ^ ( 1 + cos © ) X^ - c sin © G^ 


} 


8 


c sin © 


/ _!2 ^ \ 

\ (1 - COS ©) (1 + cos ©) J 


(2.27e) 


X Z + r X, - Y 

i 1 


{ 


( 1 - cos © ) X 


c sin © G 


} 


{ 


( 1 + cos 


© ) X + c sin © G ^ 


8 


c sin © 


{ (1 


cos ©) 


G 

? -V- \ (2.27f) 

(1 +COS ©) J 


^ (.1-003 e) _ ^ y + X Z, -r X^ 


c sin © 


( G^ Y - G Y^ ) 


(2.27g) 


(1 -I- cos ©) _ ( X Y + X Y^ ) 


-b_ 


- r X, 


c sin © 


( G Y, 


G. Y ) 


(2.27h) 


The fixed point representing rolls along the y axis are now 
represented by X,Y,Z and G as non zero and the other ooeffioients 
as zero . In this case study of KUppers Lortz instability is 
equivalent to studying the grosth of the perturbations 
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^ ^ 2 * ^ 2 ' ^ 2 ^^ 3 ' As in the earlier cases we 

again assume these perturbations to be exponentially dependent on 
time. For the onset of instability in the form of stationary 
convection ,we determine the relationship between T and e which 
makes p — 0 .The rolls along the arbitrary direction are 
represented by 

2 2 t T 

Y =(l + t^)X , G =-tX , G = r: — ^ 

2 t X 

S = b (2.28) 


Using equations (2.20a)- (2.20c) ,(2.27a), (2.27b), (2.27g) and 
(2.27h) and considering all time derivatives to be zero in these 
terms we arrive at the expression for Z^and as 





(1 - cos 0) X X, (1 +t ^ (b^^+ 4 t ^b^ ) 
bf + 4 t^ b^ - 2 r ( 1 + cos 0 ) 

**1“ “f’ 


(2.29a) 


(1 + cos 0) X X^ (1 + t^) ( b^ + 4 t^ b_ ) 
b^ + 4 t^ b - 2 r ( 1 - cos 0 ) 

(2.29b) 

The condition for the onset of Kbppers Lortz instability is given 
by considering = 0 . 
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(1 — cos 0 + c t sin 0 ) 4- ( 1 + cos 0 - c t sin 0 ) Z 

1 2 


+ ( 1 + ) c t sin 0 1^ 


2 b^( 1 + cos 0 ) 


2 b (1 - cos 0 




(2.30) 

Considering 0 = 0 in equation (2.30) results in a negative 
value of p which signifies that the perturbations coincident with 
the basic pattern decay with time and hence the original set 
renains stable to self- perturbations. The onset of Kuppers Lortz 
instability is given by putting p = 0. A plot of T against 0 
determines the domain of Kdppers Lortz instability in the T-0 
space. The minimum of such a plot corresponds to the critical 
Taylor number at which rolls along the y axis are destabilized by 
rolls at an angle of about 58°. The minimum so determined is 
approximately 23 rr* . For T greater than this value , the original 
rolls become unstable to perturbations within a certain range of 
angles with the range increasing with the rotation speed (fig. 3). 
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CHAPTER-IJT 


KOPPERS LORTZ INSTABILITY IM PRESENCE OF A MAGNETIC FIELD 
3.1 Introduction 


Rayleigh - Benard convection in a conducting fluid in the 


presence 

of magnetic field 

and 

rotation 

has 

an interesting 

feature . 

Acting independently. 

both 

rotation 

and 

magnetic field 


have an inhibiting effect on the onset of convection. When present 
simultaneously, the effect is to lower the threshold below the 
simple superposition. In this chapter we investigate the KUppers - 
Lortz instability in fluids of high Prandtl number in the presence 
of magnetic field. 

3.2 Hydrodynamic Equations 


The Navier Stokes equation has to be augmented by the 
Lorentz force arising due to the magnetic field. The effects of 
displacement current are consistently ignored. 


We have 




V + ( v»V ) V 


-♦ 

7 


P 


1 

2 


O X r p+ 


8np_ 




+2^xS+5^^(H.V,H (3.1) 

6p is the variation of density from the conduction state value, o 
is the rotation speed along the z - direction. H is the magnetic 
field in the fluid, P is the permeability and v is the kinematic 
viscosity. As in the earlier cases, Boussinesq approximation is 


assumed . 
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As the heat diffuses, the temperature T changes at various 
points . The diffusion equation for the temperature variable is 


^ + ( v‘7 )T = « V ^ T (3-2) 

with K as the thermal diffusivity. 

The equation of motion for the magnetic field is 
derived using the Maxwell's equations. Also a fluid with velocity 
V, placed in a magnetic field H, experiences an electric field 


E + ju V X H 


The current density 


Flence 


J - E + M V X H ) 
1 

E = ^- /JvxH 

t 


Using 


-> 

Curl H 


-*■ 

4 n 3 


B = - kI xW 

4 n <7^ 


we have 

Taking Curl on both sides of equation (3.5) we obtain 


Curl E = Curl [ ] - 

From Maxwell's equations we have 

-► 

Curl E = - 
-> 

■ div H = 0 

Using (3.7) in (3.6) we obtain 


M Curl ( V X H ) 


-> 

d H 


d t 


■* 
& H 
9 t 


2 ■> -> 

S. _ ^ Curl ( V X H ) 

4 n 


(3.3) 


(3.4) 

(3.5) 


(3.6) 


(3.7a) 

(3.7b) 


(3.8) 
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£JH 

3 t 


+ 


vV ) H = » 7 



( H .7 ) V 


where r? = — -i 

4 n fj <y^ 

The steady state is defined by v = 0 


(3.9) 



-> -♦ 
H = H 

o 


with as the constant magnetic field that the fluid is subjected 
to. For nathenat ical simplifications, we assume to be along 


the z direction. 

As in chapter II we scale 
i) all distances by 'd ' 
hi) time by d^/ u 

iii) temperature by AT 

iv) velocity by d/ « 

1<S 

v) magnetic field by - 


with T = , 7^^ = 

« d ^ 

and Q = 

4 n p iJ V 


the equations in the dimensionless variables reduce to 

7 ^ V 2 J w = B 0 -r +-|[ 7x7x(v7)v ] 

(3.10) 
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d t 


V 2 ] C 


T +Q ^ J ^ 

d Z ^~z ^ 




X (v’V) V 


(3.11) 


d „ 2 

<=■ 3-t - ’ 


1 

j 0 = W - ( V‘7 ) 0 


(3.12) 


^ o 2 

2 a t 




(v.7)h - (H.7)w I 


(3.13a) 


f „ ^ - V 2] . - 

3 C 


H r 

!_ ‘'j d t J ^ - 

^ z 


~\ 


— I (v.7)3 - (H.7)C [■ (3.13b) 


The strength of the non-linear terms are now discussed. In 
equations (3.1) and (3.13) the non-linear ( H*7 ) H term is of 0 ( 
~) as also the non-linearities ( v.7 ) H and ( H.7 ) v. Under all 


terrestrial conditions / a << 1 and therefore the above 


non-linear 

terms can 

be safely 

ignored . However 

the 

non-linear 

term in the 

equation 

for 0 is 

0 (1) therefore 

it 

has to be 


retained. In this chapter we consider fluids with infinite 
Prandtl number, thus the ( v.7 ) v term can also be ignored since 
it is of the 0 (<=' ”^). 

Defining J = 7 x H and for the z-components of velocity, 
vorticity, magnetic field and current density we use the notations 
w, C, h and 3 respectively. 

Assuming the time dependence of the form e , for the 
perturbations v, © and H, the conduction state (v=H=0=O ) 
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Bill be unstable if Re (p ) > o. Tor the onset of stationary 

instability p = 0 and therefore all time derivatives in (3.10) - 

(3.13) can be dropped. 

Eliminating C, e and h from (3.10) -(3.13) we obtain the 

equation 




Q ^ _ V ^ 

j* 2 

CT Z 




.« 2 f«^ _v 

^ z ^ d 


] - V 4 (-£ 


a z‘ 


j w = 0 


(3.14) 


3.3 Boundary Conditions 

If the medium adjoining the fluid is assumed 
non-conducting, then 

H = H*''' , J = 0 

z 

^ K 

where H is the field appropriate to a vacuum such that 

H = grad v# , 7 v' = 0 and on the boundary, the field is 

continuous. 

However if we consider the medium adjoining the fluid as a perfect 
conductor then no magnetic field can cross the boundary and 
therefore h = 0 

If rigid boundaries are considered then = 0 as u and v 

are ssero I and 1 = E = 0 on a boundary adjoining a perfect conductor 

K y 
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Using div J = 0 and the fact that = j =o, we conclude that 

5 J 

3.4 The onset of convection 

In the z-direction, sin nz is the lowest node that 

satisfies the idealized boundary conditions- Thus the operator 

^ 2 

^ 2 

~~~~2 replaced by — tt , Assuming the laterally infinite 

^ z 

system » the convection pattern is periodic in x~y plane with wave 

number a 

Therefore, ^ can be replaced by -a ^ and v ^ by -(n^ + a^). 
Frou equation (3.14) 


. 2. 2.3 „ 2, 2 2. „ 2 

R(a> = - 1 + 

a a a 


Q T 

~2 2 2 ^ 
a [(7T^+a^)+Q77^] 


(3.15) 


Thus we see that the stabilization by Q alone is 
expressed in the second term and that stabilization by T alone in 
the third term. Lowering of the Rayleigh number below the 
threshold is exhibited by the fourth term. 

Minimizing R(a) with respect to a gives the critical wave number. 
The equation that 'a^' satisfies becomes 


3(rr^+a (n^+a Q 

c c 




7 2 2 

2a ^(^^+a + 

c c 


Q rr^ + (n^+a^^)^j ^(n 


+ Q 


(3.16) 
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3.5 Kuppers-Lortz instability 


For fluids with infinite Prandtl number, ct ^ >> a . The 

time dependence is to be retained therefore only in equation 
(3.12). 


In equation (3.13) if we operate 


a 


d 2 


on both sides we obtain. 


„ 2 a h a 
-5-i = —a" 

d z 


(3.17) 




operating on equation (3.11) gives 


-7 - T w + Q J 


a 2 


a 2 


(3.18) 


operating ^ ' on equation (3,10) 


„ 2 f a_ „ 2 
^ 1 a t 


v^*] » 


2^ 2 g 

i 


2 a c 


+ Q 7 ‘*(7 


2.„2 a h 


a 2 


a z 


) 


(3.19) 


Eliminating h by using (3.17), equation (3.19) reduces to 


7 2 1" ^ ^ 7 2 _ ^ 4j ^ 2^^2 Q _ 


2r, 2 ^ _r7 2 a C „ „2 d w 


- Q 7" 


a 2 


a 2 


(3.20) 


Consider 7 of equation (3.18) 


7 2w + Q -t— 7 2 j 


-7 ’ C = 


a z 


a z 


(3.21a) 


Using equation (3.14) 


2j _ J_L_ equation (3.21) becomes 
y tM a z 


-7 ^ c 


a 


a 2 


^ ^ 

7 ^ w - Q ^ 2 

a z 


(3.21b) 
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or. 



C + T 


d 


d 2 


V" w = 0 




) 0 + w 


-> -> 

(v*v ) e 


(3.22) 

(3.23) 


Thus the systen has been reduced to three variables. 
Expressing w, C, 0 in vectorial form. 


X = 


C 

0 


a compact form for expressing (3.21) - (3.23) would be 


X + (AR) L X = N( X, X') + O' I- X 
* . a t I 


13 


(3.24) 


N( X, X ') is the column vector 


0 
0 

-> 

^ (v *7)0' 


AR = R » R 

€ 

is the operator 


V V ■*- Q D^) 

T D V ^ 

1 


-T V ^ D 

V Q 
0 



(3.25a) 


and L is 

% 


r 0 0 7 ^7 2 1 

i 

0 0 0 

, 0 0 0 . 


(3.25b) 


Assuming that X represents the cylindrical rolls that 
are formed at the onset of convection, using perturbation theory 
by a power expansion in AH , it is possible to determine this 
state for slightly super - critical R . 
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Expressing the deviation in R from 


R as AR we 


00 


AR = 2 

n=l 




R 


n 


-> 00 ^ 

X = E X 

n=l 


« represents the anplitude of the non-linear convection 
zeroth order in e 




> H(X'.X) = 

' 0 1 
0 


L 0 ^ 

A 


0 J 


AR 


R 


Equation (3.25) becomes 


^-0*0= - K ^ , 


We obtain at the zeroth order 


= cos a X sin nz 


u. 


n 


sin a X cos nz 


2. 

(n +a ) 
nr '■ “ 

a ,2 2^2 2 „ 

(n +a ) +n^Q 


sin a X cos nz 


TT T (n^+a 

^vr2^„ 2.2^ 2 q 

(rr +a ) +n Q 


cos ax cos nz 


0 


(rr^+a 


sin nz cos a x 


expand 

(3.26) 

(3.27) 

At the 

(3.28) 

(3.29) 

(3.30a) 

(3.30b) 

(3.30c) 

(3.30d) 

(3.30e) 


At the first order equation (3.25) becomes equivalent to 
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1 


' 0 0 ■ V 2 . 

1 




0 

1*0 



0 0 0 



— 

0 


1^, J 


.0 0 0 ^ 


o 

e 

^ o 


-> 

, (v *7)0 ^ 


Applying the solvability criteria described in Appendix B, we have 

^ ^ol ^ ^ ^ol ^ = 0 (3.32) 

The first expression gives a zero contribution while 
1 2 

f not equal to zero. This necessarily inplies that 

= 0 . We also have w^= = u^= v^ = 0 (3.33a) 


sin 2^z 

o y 2 2. 

8 (rr + a ) 


At the next order , 
X, + R L X, + R 

0 2 i 1 i. 


the equation becomes 
^ Xo = N (X^,Y^) + H (X^.Y^) 


(3.33b) 


(3.34) 


Using equation (3.34) ,we arrive at the solvability criteria as 
- < X^l > + < X^l N (X^, Y^) > = 0 (3.35) 


This results in R = ^ « — 2 

8 ( n'^+ a ) 

Since R,> 0 ,it implies that the bifurcation to the convection 
state is positive for all values of t and .0. Kdppers and Loft 
have found the bifurcation to be greater than zero for Q = 0. This 
result is thus unaffected by the presence of the magnetic field. 

The investigation of the stability of the cylindrical rolls 
against perturbations by another set is the basic feature to be 
explored for the study of Kdppers Lortz instability. We assume a 
perturbation vector Y with time dependence of the form e .A study 
of vrowth of this oerturbation determines the stability of the 
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original set described by X. With operators and l^as described 
above , the hydrodynanic equation (3.21) - (3.23) can be 

represented as 

- p a H Y + L^Y + ARLJ = N(X.Y) + H(Y,X) (3.37) 

where the matrix M is 

r 0 0 0 . 

0 0 0 

^0 0 1 J 

He expand Y exactly as we represented X .Thus we have 

to CO 

' ’'n -1 P = ^ / P„ <3-38) 

n = 1 n = 0 

As we are considering the stability of the marginal state^, at the 
zeroth order we choose p^ = 0 .If the perturbing roll is also a 
cylinder with the axis along the y direction ,the original roll 
remains stable. 

At the zeroth order we have = 0 


{ 


w = Re < e 1 + “r’'’ Sin az 

0 


} 


(3.39a) 


u,. = Re 


[{ 


irerk^C^ + a ; 

( (n^+ ^ Q 2 


ink. 


} 


i (k X +k_y) 


e 


1 .-^ “ 2 '' cos nz 


(3.39b) 
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V, = Re 
o 


r / i ^ T (n^ + a^) 

i n , 

^ ^ r /rrZj 2 .2 2 2 

E (n + a )^ + Q ] 

2 1 

a 

i (k X + k_y) 1 

e * ^ J 

cos nz 


(3.39c) 


<9, = Re 
o 


i( k^x +k2y) 


^2 ^ 2 
Tx + a 


sin nz 


(3.39d) 


= Re 


^ « i (k X + k y) 

n T ( + a^) e 


(n 


2 


2 2 2 
+ &r + Q 


cos nz (3.39e) 


At the first order in e , the equation reduces to 


= ^ <^o'Yo) + 


(3.40) 


Solvability condition becomes equivalent to 

< Y^l M (X^,y^) + H (Y^,X^) > + < Y^l c p y ^> = 0 (3.41) 

In equation (3.41), since the first tern gives a zero 
contribution and since <Y^| M Y^> is not zero it necessarily 
implies that p^ = 0. 

Thus to investigate the Xdppers Lortz instability, terms of 
higher order in « need to be considered. With p^= 0 ,R^= 0 we 

obtain expressions for w^,u^,v^,C^and as 
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[{^ 


"t = I 1 *+ ® 


iC k^+a )x 


i( k^- a )x ^ ik,y 
e 


H- i e ‘ - j ] 


Sin 2 nz (3.42a) 


C. = Re [ [b 


i(k + a )x 


i(k^- a )x >1 ik^y 
e 


+ B e ‘ 


COS 2nz (3.42b) 


= Re[ ( a 


i(k + a)x 


i(k -a)x >, ik y 
+ C e " \ e ^ 


} ^ " ] 


sin 2nz (3.42c) 


u, = He 


{(■ 


(2ni (k^+ a ) + ik^ ) 


(k^+ a )^ + k/ 


i(k +a )x 
1 


j. (2rTi (k^- a ) A_ + ik^ B_ ) ^ 

\ ^ J 


(k^- a )^ + k/ 


ik^y 

e cos 2nz 


i(k - a )x 
1 


} 


(3.42d) 
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= Re 


[([ 


2ni - i(k^+ a ) 

i (k^+ a )x 

2 2 
(k + a )'^ + k^ 

1 2 

1 6 

2n:i k_ A - i(k - a ) B 

r 

1 - 7=^ -1 

^ .2 .2 

J ® / 


i k-y 1 

e cos 2nz (3.42e) 


where A. 


2 (1 ± cos 0) R n (1 + cos ©) 

(A a^) 


(3.43a) 


■2n^R T (1+cos 0)[2a ^(l±cos e)][2a ^(l±cos ©)4-4n ] 


± “ 2 , 2 

n +a 


2 . 2 .. 2 , 


C{2a ^(l±cos ©)+4n } +471 Q] A+ 


(3.43b) 


»n(l+co 3 e)f {2a ^(l ±cos e)+ 4 n^}^+ 4 rT^Q{ 2 a ^(l±cos e)+4n }] 

" (n\a ^)[2a ^(l±cos ©)+47 t^] A^ (3.43c) 

^ ^ ^ ^ . a _2 T 2 


((t+cos 0) 4n^T [2a ^(l±cos ©)+4n ] 


2 . .._ 2 , 


(rr^^a ^)C2a ^(Itcos 6 »)+ 4 n 23 A, [{2a (l±oos a)+4n }+4n Q] 


1 2 a ^<l±cos 0)+4n^3^ +4n^Q[2a (l±cos ©)+4n ] 

4n^C2a ^(l±cos e)+4nh^ _ 2a 2 r^(1±cos ©) 

C 2 a ^(Itcos ®) + 4 n^ 3 ^+ 4 n Q 

* —j—.. cjA have froBi Bq. (3.37) 

At the second order , we 


(3.43d) 
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(3.44) 


and applying the solvability criterion arrive at 


, sin & 

B 

+ 


1 uTsin 0 (n^+a ^)[C^-C_3 

t 4(u^4a ^) 

- 2( 1+cos 0) 

2( 1-cos 0) 

.«r/n2^ 2.2,„2m 

4[(n +a ) +u QJ 


+ 3 - [(1+cos 0 ) C_ + (1-cos 0) C 




f 3.451 


For T =0, P 2 is negative as it should be and increases as t 
increases, having a zero, when 


sin 0 r 



] + ”[ C^d-cos 0) + C_( 1+cos 0)j 

^2.^ 2 [ 2( 14 CO 

s 0) 

2(l-cos0) 

2 7 

nr (n +a ') 

sin 0 

r 1 

(3.46) 

+ - 2 2 2 

(u' »a 

jrr^« 

■ [ ] 

= 0 


The condition for Rdppers-Lortz instability is thus given by 

Eq.(3.46). 


3.6 KOppcrs I. ortz instability with Galerkin trunca tion, 

»e now study the Rdppers Lortz instability using a Galerkin 
model. The minimal model is considered first for the calculation 
of the critical rotation speed. A more complete representation is 
then considered which though complicated, has the advantage 
yielding realistic answers. For the stress free boundaries we 
consider the following representations for the rolls along th y 

axis 

w = a<t> cos ax sin nz 


C = f<t) cos ax cos nz 


(3.47a) 

(3.47b) 
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& = b(t) cos ax sin nz + c(t) sin 2nz 

h = h(t) cos ax cos nz 

3 = a(t) cos ax sin nz 


(3.47c) 

C3.47d) 

C3.47e) 


Introduction of this model in the hydrodynamic equations results 
in a set of 6 ordinary differential equations corresponding to the 
time evolution of the coefficients of the Fourier modes. With 
X,Y,Z,G,H,J as the scaled version of the coefficients 
a( t ) , b( t ) , c( t ) , f ( t ) , h( t ) , j ( t ) respectively, we have the Lorenz 
model as 


X = o- 


{ 


X + Y + tG - 


2 

Q 


/ 2 ^ 2.2 
(n + a ) 


« } 


<3. 48a) 


Y + r X - X Z 


(3.48b) 


= X Y 


b Z 


(3.48c) 


{ 


= - o- ■{ G + tX + 


Q J 


2.2 

(rr + a ) 


} 


(3.48d) 


c/ H = X - H 
2 


(3.48e) 




(3.48f ) 


The stationary state solution of these equations is represented 
by setting all time derivatives to be zero. This gives X = Y = Z = 
H = J = G = 0. This represents the state of no convection . However , 
this equilibrium state of the fluid is perturbed at some critical 
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tenperature gradient which is determined by the linear stability 
analysis about the conduction state .As in the earlier cases , we 
superimpose x, y, z, g,h and j as the perturbations, on the 
conduction state and set up the linearized hydrodynamic equations 
in these perturbations. Assuming time dependence of the form of 
e^*^ , the marginal state will be defined by Re (p) = 0. In general, 
the marginal state can be either stationary or oscillatory 
depending on whether Im (p) = 0 or not for Re (p) = 0. For the 

conducting fluid placed in a magnetic field the marginal state 
can be oscillatory only when the magnetic Prandtl number 
defined as is larger than the thermal Prandtl number ( ct 

This condition is equivalent to having > 1 - For 

->3 

terrestrial conditions , this is of the order of 10 . Hence we 

expect that the onset of instability is via a stationary state for 
fluids under consideration. Thus the critical Rayleigh number for 
the onset of convection is determined under the condition Im (p) = 


0forRe(p)=0 . 

The disturbances are thus temporarily guided by 


- 1 


1 

0 


0 -1 


0 


0 




0 


, 2 ^ 2 .2 

(n + a ) 


0 


0 


0 


tt^Q 


(."4- a^ )2 


-1 


0 


0 


-1 


= 0 


I: 






I 

2 ! 
yj ! 
u ’ 


(3.49) 
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With the assunption of non zero values of the variables x,y,g,h 
and j . We notice that the equation for Z decouples when linearized 
in these perturbations. Equating the determinant of the above 
matrix to zero determines the critical Rayleigh number for the 
onset of stationary convection. 


R 


, 2 2.3 


2 2 2 
(n^+ a^) 




Q T n 




(3.50) 


with r 


R a" 


, 2 2 .3 

( fi + a ) 


(3.51) 


The second and the third terms represent the inhibiting effect of 

the magnetic field and the rotation respectively. The last term 

expresses the reduction in the Rayleigh number when the fluid is 

subjected simultaneously to the two constraints. The critical 

wave number a is determined by minimizing R(a) with respect to 
c 

a. Thus a satisfies the equation 


3 ( + a^ )^ 

C 


, 2 . 2 .3 

( n + a ) 


Q n 


T n 


T Q 


4 r 2 a^ (n^+ a^) + Q + ( rr^ + a^ )^ 1 

n C 


0 

(3.52) 


The critical Rayleigh number and the critical wave number 
j - 4- r> are in complete agreement with the results arrived at 



by considering perturbation techniques. Thus we conclude that our 
truncated systen is an accurate nodel for the investigation of 
Kuppers Lortz instability . 

As the critical Rayleigh number depends only on the magnitude 
of a , it is obvious that an infinite degeneracy exists. All steady 
state solutions which have the same wave number but different 
directions are equally preferred due to the isotropy in the 
horizontal directions. In order to consider which of the solutions 
will be physically realized ,we need to enter the non-linear 
regime. The steady state rolls that are formed at the onset are 
described as 


Y = 


1 + 


1 + 


Q 


. 2 ^ 2.2 

(n + a ) 


7T^ Q 


,2 2 ^ 2 

+ a ) 


} 


(3.53a) 


G = 


- X t 


1 + 


2 

Q 


, 2 ^ 2.2 
( rr + a ) 


(3.53b) 


H = X 


(3.53c) 


J = G 


(3,53d) 


The steady state rolls assumed along the y axis which are 
described by equation (3.53a) - (3,53d) are then perturbed by 
rolls along an arbitrary direction. Thus we need to include the y 
dependence as well in the expansions for w, C» h and j. We 



therefore represent these variables as 


- a(t) cos a X sin nz + a^<t) cos ( k^x +k^y ) sin nz 


C - cos nz f(t) cos a X + f^(t) cos (k^x +k^y) ^ 

<9 = sin nz | b(t) cos ax + b^(t) cos (k^x + k^y) 

+ sin 2nz | c(t) + c^(t) cos (k^x +k^y + a x) + 

+ c^Ct) cos (k^x +k 2 y - a x) I 


h = cos nz 


j = sin nz 


h(t) cos a X 


I j(t) COS a X 


+ 


+ 


h (t) cos (k X + k y) 

X ± Z 


j^(t) cos (k^x + k^y) 


} 

} 


Introduction of these expansions into the equations (3.21) 
results in a 13 x 13 Lorenz systen. 

The rolls along an arbitrary direction are described 


{ 


1 + 


n^ Q 


1 + 


n^ Q 

(n2 . a2)2 


/„2 ^ 2 
(n + a 


G. 



t 

Q 


(3.54a) 

(3.54b) 

(3.54c) 

(3.54d) 

(3.54e) 

- (3.23) 

by the 

(3.55a) 


(3.55b) 
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(3.55c) 


J = G 

1 1 (3.55d) 

since these rolls satisfy exactly the sane equations as the 
original rolls. The condition for Kdppers Lortz instability for 
the nininal nodel is given by 


t sin ® - 


b^( 


b^(l+ cos 0)^+b (1-cos 0)^ 


1+cos 0) - b (1-eos 0) 


(3.56) 


where b^ and b_ are defined as b^ = 


4n^ + 2a^ (1 ± cos 0) 


2 2 
TT + a 


(3.57) 


For Q = 0 , this equation becomes identical to equation (2.26) From 
this expression we see that T is a monotonically increasing 
function of Q . The critical 0 for the onset of Kdppers Lortz 
instability is found to be about 58°. It is very weakly dependent 
on the magnetic field . 

The minimal model that has been discussed above has the 
advantage of being handled easily . It is able to show all the 
qualitative features of Kuppers Lortz instability in the presence 
of the magnetic field. However to obtain more accurate results we 
need to consider a higher order truncated model. We introduce 
additional terms in the w , C, h, j expansions. Accordingly we 


have 
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w sin nz a(t) cos a x + a^(t) cos (k^x +k 2 y) + sin 2nz 

cos (k^x + a X +k 2 y) + a^Ct) cos (k^x - a x + k^y) 


(3.58a) 


C - cos nz f(t) cos ax + f^(t) cos (k^x + ^ 2 ^^^ 

jf^Ct) cos (k^x + a X + k^y) + f^<t} cos (k^x - a x + k^y) | 


(3.58b) 


® - sin nz ■|^b(t) cos a x + b^(t) cos (k^x +^ 2 ^^^ ^ 2nz 

|c(t) + c^(t) cos(k^x + a X +k 2 y) + c^it} cos(k^x-a x + k^y)! 


(3.58c) 


h = cos nz l^h(t) cos a x + h^(t) cos (k^x +^ 2 ^^^ ■*■ 

|h^(t) cos (k^x + k^y +a x ) + hg(t) cos (k^x +^ 2 ^“ ^ 

(3.58d) 

0 = sin nz -|o(t) cos a x + j^(t) cos (k^x +k 2 y)|’ + sin 2 rT 2 

j 2 (t) cos ( + k^y +a x) + j 3 (t) cos (k^x-a x+k^y)!^ 

(3.58e) 


Substitution of this nodel in the hydrodynamic equations 
results in a dynamical system of 21 equations. 


= <y i -X + G 


t + T - . « 1 

^ 2 .2 J 

(n + a ) 


Q H 


r u « n 

] -X^ + G^t + Y^- 2-2- \ 

I ^ ‘ ( 77 ^ + a^ )^ J 


(3.59a) 


(3.59b) 
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K 


2 t G, 2 rt"* Q H 

X. = ^ 1-. X, . — ^ _ 


+ a2 )2 J 


(3.59c) 


{-bj 


2 t G. 


= O' i- b_X3 + ^ 


2 rr^ Q H 

■ 7 ^ 


Q H 3 . 

. a2)2 ; 


(3.59d) 


G = 


G. 


a 


I - 6 . tX - I 

^ (n^+ & J 


r a 1 


(3.59e) 


(3.59f ) 


{ 


2 n'- Q , 

+ ^ <.2^ a" J 


O' -i 2 t X^ + Sj b 


(3.58g) 


{ 


O' 2 t X, + G^ b + 

3 3 — 


« J 3 \ 

^ a2 )2 / 


2 Q J 
(n^ 


(3.59h) 


Ho = X - H 
2 


(3.59i) 


H o = X - 

* 1 211 


(3.59j) 


H, -a = - »+ ». + 2 X, 


(3.59k) 


H 0 = - b H„ + 2X 

**3 2 — 3 a 


(3.591) 
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J = G ~ j 


(3.5911) 


J a = G - J 
12 1 *^1 


(3.59n) 


J<y=2G-bJ 

3 2 2 - 3 


(3.59p) 


Y = - XZ+rX-Y- 




( 1- cos 0 ) + cG sin d ) 


( 1+ cos 0 ) - c G^sin B ) 




c sin & 


{ 8(1 


cos O ) 


I 

( 1 + cos 0 ) J 


(3..59q) 


Z + r X, - - 




cos 0)X - c G sin 0 


(1+ cos 0)X + c G sin BY + Y c sin 0 


/ 

18(1 


cos B ) 


8(1+ cos 0 


(3.59r) 


Z = XY + X^Y^- bZ 


(3.59s) 


( 1 - cos B ) 
2 


(X Y, + X Y) - b Z, - r X, + 

X 1 *1" 1 ^ 


c sin 0 


( G^ Y - G Y^ ) (3.59t) 
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Zj, = — ( 1 + cos 0 ) ( X Y) - b_ - r Xg + 

— - - — ( G Y^ - G^Y ) (3.59u) 

We have rescaled time by a factor of o’. 

Since we require to investigate the instability right at the 
threshold we drop all the terms containing the time derivatives. 
The condition for Kdppers Lortz instability is then obtained by 
substituting the values of the variables from equation (3.59) 
after equating the time derivatives to zero in the expression for 
Y^ in equation (3.63g). 


2 sin^0 R a^ 


2 sin^0 R a^ 


With X, = - 
2 


and Xg = - 
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value of about 60 The results that are obtained using the 
Galerkin truncation are in agreement with that obtained in section 
3. 5 using the perturbation techniques. 

3.7 Kuppers Lortz instability in finite Prandtl number fluids : 


In this section, we lift the constraint of o' >> 1. For fluids 
with low Prandtl number, the onset of convection is oscillatory 
when the fluid is rotated at a sufficiently high speed. However, 
instead of rotation, if we subject such a fluid to a magnetic 
field, the condition for the onset of convection as oscillatory is 
that the thermal Prandtl number has to be small in comparison to 
the magnetic Prandtl number. However under all terrestrial 
conditions, the ratio of / o' ) is 0( 10 and thus the onset 

of convection as stationary pattern is ensured. To investigate the 
Kuppers-Lortz instability in fluids with low Prandtl number, we 
subject it simultaneously to rotation and magnetic field. The mode 
of instability at the onset is now a complicated function of the 
Taylor number, T, the Chandrashekhar number Q, thermal Prandtl 
number o' and the magnetic Prandtl number o'^ . We remain in that 
region of Q where the instability as stationary convection is 
ensured and in particular we investigate Kuppers-Lortz instability 
in mercury. 

The Hydrodynamic Equations 

The non-linear terms in w and C can no longer be dropped 
since they are of 0( o' ^ ) . As the non-linear terms in the 
equations for h and C are of 0( ^ ) and in mercury ( o'^ / ) 

is of the order of 10 these terms can be neglected. The 
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non-linear term in the equation for & being 0( 1 ) is retained. 

The governing equations for the finite Prandtl number becomes 




1 

O' 


7 X 7 X ( vo7 




T 


c 

d z 


+ 


Q V 


2 ^ h 

d z 


(3.60a) 


( K H ’ 
[ - K - ’ " ] ® = " 


X ( V*v 


- ( vv 


d w 
d z 


^ c 

d z 



T 


^ W 

^ z 


+ Q 


^ j 


(3.60b) 

(3.60c) 

(3.60d) 

(3.60e) 


The subscript z denotes the z-component of the additional 
non-linear terms. 

We study the KOppers Lortz instability via the Galerkin 
truncated model used in Section 3.6. Introducing these 
representations in the governing equations result in the following 
system of equations. 
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X = cr 


f-X + Y+ tG ^ ^ y H 1 - I r — J 

(n^+a 2)2 J 4 L („2 


^ 2 
+a ) 


cX.sin O 


{- 


c (l~COS 0) „ (1+cos Q) 

2 (1+cos 0) 3 (1-cos 0) 




G (1-COS 0) +G ( 1+cos 0) 


} 


-X X ( 1 - cos 0 ) - X X ( 1+cos 0 ) - c sin 0 G^ ( X - X ) 

* ^ jL 3 X 2 3 


-X 


sin 0 


{ 


2 


G 


i 2 1 (1+cos 9) (1-cos 9) 


}] 


(3.61a) 


X. = - 


{ -X. + Y, + t G, - -27^ - S 

V (tt +a ) ^ 


n 


2 2 
(^ +a 


cX sin 0 


{- 


p (1-cos 0) p (1+cos 0) 
^^2 (1+cos 0) *^3 (1-COS 0) 


-J+ G G 2 ( 1-cos 0) + Gg( 1+COS 0) ^ 


-X X,(l-cos 0) - X X,(l+cos 0) + c sin 0 G (X,- X,) 

2 3 ^3 


+ X 


sin 0 


{ 


(1+cos 0) (1-cos 0) 


} 


(3.61b) 


X = cr 
2 


{ 


^ ^ ^2 2 Q 


2 z 


H 


,2^^ 2.2 “2 b. 
(rr +a ) + 


(1+cos 0) 


} 


4 b 


2 2 
8 ^ ^ _ ^x2 . 8 n 

2, X ,_2._ 2. 


c sin 0 (1-cos 0) 


(" +a *') 


(n^'+a ‘') 


(G X^ - G^ X) - 8 sin^ 0 G G^ + b^ | 4 X X^( 1-cos 0) 


+ 2c sin 0 (G^ X - G X^) I 


(3.61c) 
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r ^ ^ % 9 a 2 z, . 

*5 + -bT- ' 2j2 "a - SZ- } 

1 r ® o 2 

4~b ~ 2 ^ ^ ^1 (1+cos 9) - — ^ ^ c sin 9 (1+cos 9) 

~ (n +a ) (n^+a ) 


(G - G^ X) - 8 sin 0 G G 




4 X X^( 1-cos 9) 


-2c sin 0 (G X - G X ) 
1 1 


(3.61d) 


a = cr 


{- 


t X - 


rr^ Q J 
(rr^+a ^ 




G^d-cos 0)+ X^GgCl+cos 9) 


. c sin 0 ^ c sin 0 G^G_(l-2cos 0) 

2(l+COS 1 2' ^ 2(l-cos O) 


+ — 2 ^ ^ ^ 

(rr +a ) 




Y (1-cos 0) Y Y (1+cos 0) 
^^2 (1+cos 0) 1-^3 (1-cos 0) 


+ GJ X^d-cos 0) + Xgd+cos 0) 


(3.61e) 


G, = o- 


{ -G, - t X. 


2 

rr*^ Q 
2 

(JT +a 


)2 


X G^d-cos 0) 


+ X G^d+cos 0) - 2jho^l h { ® 


c sin 0 ^ ,1 o ON ^ o / V V d-cos 0) 

■*■ 2( 1-cos 0) ^ ^ ( 7 T^+a \ ^ ^2 (1+cos 0) 


* *3 ( i^cOS I ) } ®> + X3(H-OOS 9) J 


(3.61f ) 


G = <y 
2 


2 

I -G^ - 2 t X^ ^ - ^ } (3-61g) 

V (n +a ) 
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G = a 

3 


{ -s, b_ - 2 t X, - - y! - a J I 

(n^+a ^ ) 


(3.61h) 


H 

2 1 


- H. 


(3.61i) 


H 


X - H 


(3.61j) 


Ha 


2 3 


J = 


2 H^ 


2 X 3 - b_ H 3 


G - J 


(3.61k) 

(3.611) 

(3.61in) 


O' J = G - J 

2 i 11 


(3.61n) 


O' J = 
2 2 


2 G^ - J, 


(3.6I0) 


o- J = 2 G - J 

2 3 3 3 


(3.61p) 


Y = -X2 + rX-Y-X Z ^ ^ - X Z 

11 4 12 4 


^ « Y c sin © 

^ ^ ^ p / 7 _ 7 ^ 

4 112 '' 4 1 2( 1-cos ©) 2(l+cos ©) 


{ 


} 


(3.61q) 


Z = XY + XY-bZ 

1 3 . 


(3.61r) 


X Z 


X Z 


Y = - X Z + r X - Y^ 

t X , 1 1 


(1-cos ©) - 


( 1+cos ©) 


O' G" 

. c sin , ^.Yc sin 0 f 3 2 

+ 4 It - -62 ; + 4 2 (l-cos ©) 2 ( 1 -cos ©) 


} 


(3 .61s) 
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2 . = 


(1 . . 90 S e) ( Y X + X Y > - b. Z + ^ (GY -GY) 

^ j- 1 +1 Z i 1 


- r X 


(3.61t) 




(l+cos . ( Y X^ + X Y^) - b_ - — - -|- --- iGJ - G Y^) 


r X. 


(3.61u) 


In the above equations we have used 




-> -+ 

V X ( v-v ) 


a . 3 V. 3 V. 

( 5-^ yX } - 

Jz ^ J 


■* -> 


V ^ (V*^ ) (3.62a) 


and 




-*■ 

X (v*'7 ) V 


I - (^.<7 ) ( - (W7 ) w 


(3.62b) 


Substituting for Y, Z, G in the expressions for Z^and Z^ we 


arrive at 


X X 

Z^ = (1-cos 0) 




1 + 


1 +n^Q 


t^ , Q _ \ 

( n 2 ^a 2)-2 ^>2 / 


(3.63a) 


Z^ = ( 1 +cos 6 ) ^ 


b- I 


1 + 


1+a^Q (7i^+a ^) ^ ( 


rt^Q \ 


r X 3 
b_ 

(3.63b) 


FroB equations (3.61) - (3.63) we have 
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If. = - z, { 


1 - cos O + 


l+rr 


c t sin 6 1 rr r ^ 

2q <.2+a 2)-2 1 ■ I 


0 


c t sin © 


2 2 2 
1+rr^Q (7T^+a ^ 


)-2 } 


+ ct sin © 


/ !i!a_ 

I 2) 


2 2 2 2 
l+n‘^Q(n^+a 


7 ^} 


b, X 
+ 2 


b_ X. 


. 2 „2r,,„2. 2>-2 (1+cos ©) 

b^ +4 n Q(77 +a ) 


2 2 2 2 
b_ +47T^Q(77^+a 


I 

^-2 (1-cos ©)J 


Where X^ and X^ are given by 


(3.64) 


N. 


O' 


X_ = 


^3 = 


M 


N. 


a 


D_ 


M = “2 sin B i 1 + 


{ 


Q 

(n2 . a2)2 


1 + 


tt^Q 


2 2 2 
(n^+ a^)^ 


} 


R = 

2 


2 n^(l - cos ©)^ 

(7T^+ a^) 


- b^(l 


cos ©) - 


2 2 
2 sin'^ © 


fi - / U] 

^ (n + a ) 


D. 


b? + 


2 

4 t/ b, 


4 n^Q 


4 Q 

. 2^ 2.2 
(n + a ) 


2 r (1 + cos ©) 
b. 


(n‘'+ a^*) 


2:2 + 
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N = 

3 


2 (1 + cos 6>)^ 

( TT^ + ) 


- b (1 + cos 0) - 


2 sin^ 0 
Q 


[ - T ^ 


D 3 . b: 


2 

4 t b 


2 

4 rr^Q 


/ 2 ^ 2.2 
(rr + a ) 


+ b" 


4 7T^ Q 

(.2 ^ a 2)2 


2 r (1-cos &) 
b 


The condition for the onset of KUppers Lortz instability is 
given by Y^= 0. The effect of including the non linear terns in 
the equations for the velocity and vorticity is to reduce the 
critical Rayleigh number ( fig. 6 ). The critical angle is weakly 
dependent on the Prandtl number . Thus the presence of the 
magnetic field has enabled us to investigate KUppers Lortz 
instability in fluids with arbitrary Prandtl numbers . 
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CHAPTER IV 

KOPPERS LQRTZ INSTABILITY IM FLUIDS WITH INSULATING BOUHDARIES 

4 . 1 Introduction 

He consider a horizontal layer of fluid rotating about a 
vertical axis to be confined between stress free boundaries. The 
thermal conditions on the boundaries are that of constant heat 
flux. We study the mode of instability that is preferred at the 
onset using linear stability analysis. We find that the onset of 
instability is in the form of long wavelength rolls, below a 
critical Taylor number. Next we address the question of Kdppers 
Lortz instability in these long wavelength rolls. Investigation of 
Ktippers Lortz instability is carried out via a Galerkin truncated 
model . 

4.2 Hydrodynamic Equations in a dimensionless form 

From Appendix A, the hydrodynamic equations in the dimensionless 
form are represented as 

V ^ ( V 2 _ ) w = T — - R ® (4.1a) 


( ^ 


2 


d 

& t 


) C 


T 


& W 

9 z 


(4.1b) 


( ^ 


d 


) 0 


( V . V ) © 


2 


-<y 


d t 


w + 


(4.1c) 
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4.3 Boundary Conditions 


We consider the fluid to be confined within stress free 
boundaries and accordingly we have 

d ^ w 

w = - 0 (4.2a) 

d zT 


■ Q C 
d z 


0 


(4.2b) 


The thermal conditions on the boundaries are that of constant heat 
flux, as we are assuming the fluid to be bounded by insulating 

surfaces. Thus unlike the case of perfectly conducting boundaries, 

where we had O (the temperature fluctuations from the conduction 

S 

state ) equal to zero, we now have ^ — =0 (4.2c) 

on the boundaries. 

4 . 4 Tli e_ LQrenz._aQdeI 

A minimal truncated model which is consistent with the 
boundary conditions is assumed to be the solution of the equations 
(4.1a) - (4.1c). Introduction of the minimal model into these 
equations results in a system of 4 ordinary differential 
equations. Irrespective of the thermal conditions on the 
boundaries , the motion of the convecting fluid is described in 
terms of the variables w , C and S . 


w = a(t) cos a x sin nz 


C = f(t) cos a X cos TTz 


9 - b(t) cos ax + c(t) cos nz 


(4.3a) 

(4.3b) 

(4.3c) 


Substitution of these 


representations in equations 


(4.1a)- 
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(4.1c) results in the dynamical system (after scaling time by <y) 


X = ^ (-X + Y + tG) 


(4.4a) 


y 


- X Z + 


8 r X 



a2 Y 


.2 + 


.(4.4b) 


Z 


X Y - 


rr^ 7 

n Z 


.2 


(4.4c) 


G = - o'(tX + G) (4.4d) 

The equilibrium state of the system is defined by the solution of 
the equations (4.4a) -.(4.4b) under the condition of all time 
derivatives being zero. The steady state is thus defined by X,Y,Z 
and G are zero. This corresponds to a state of no convection. 
However as in the earlier cases the fluid remains stationary only 
below a critical temperature gradient characterized by the 
Rayleigh number. Above R^, the conduction state gets destabilized 
in favour of a convection state. Whether the instability takes 
place via a stationary state or an oscillatory state is determined 
by the linear stability analysis. Employing this technique, we 
calculate the critical Rayleigh number for the onset of stationary 
as well as oscillatory instability. The marginal state would be 
the one corresponding to the lower Rayleigh number. In the limit of 
large o' , the onset is stationary while for o tending towards zero 
the onset is oscillatory. 
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^ ^ Linear Stability Analysis 

He perturb the conduction state (X=Y=Z=G=0) 
infinitesimally. These disturbances are assumed to have a time 
dependence of the form e^^. He set up the hydrodynamic equations 
in these perturbations (neglecting all second and higher order 
terms in them ) and determine the critical Rayleigh number and the 
critical wave number . If the onset is assumed to be stationary, 
we have p = 0 and the Rayleigh number determined under these 
conditions defines the threshold for stationary convection . If we 
assume the onset as oscillatory, we consider p = ± iw with ^ as 
the frequency of oscillation and determine the critical Rayleigh 
number corresponding to it . As we are considering the marginal 
state the Re (p) = 0 

He assume perturbations of the form x,y,z,g to be 
superimposed on the conduction state X=Y=Z=G=0 and study 
the growth of the disturbance. Linearizing in the perturbations, 
the equation for z decouples and thus only the variables x,y,g 
need to be considered. These linearized hydrodynamic equations are 


px = O’ (-x + y+ tg) 


(4 . 5a) 


P g 


- O’ (_ tx + g) 


(4.5b) 


p y = 


8 r X 


a^ y 


n 


2 , 2 
n + a 


(4.5c) 


Equations (4.5a) - (4.5c) can be represented in a matrix form 
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as 


P + O' 

8 r 

" n2 

a t 


- a 


P + a^( 


0 


2.n2 + 


- o t 

0 

p + o 


y 

g 


= 0 


For a non trivial solution of these equations, the 
solvability condition is that the determinant of the ma 
multiplying the variables x,y,g must necessarily be 
results in a cubic equation in p . 


+ p2 ( 2 <?• + 


a 


+ a^ 


) + P ( ^ 


(1 + t^) - 


2 a 


rr^ ^ 2 

n + a 


) + 


2 2 ^ .2. 

a a ( 1 + t ) 

2~7 2 
TT + a 


O’ 


2 8 r 


8 r O' 


n 


= 0 


rr 


(4.6) 


i- - (A RV is used for 
For the onset to be stationary equation 

determining the critical Rayleigh number with the condition that p 
= 0. This is equivalent to having the constant term m equation 
(4.6) as zero. The critical Rayleigh number r__ is determined to 


r = 
c 


(1 + t^) 
8 ( + a^ ) 


(4.7) 


Substituting 


^c = 


R 


(rr^ + a^ )^ 
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2 

T 


( 4- 


We obtain 


F^(a) 




2+ a2)2 + 


T 

a^ 


} 


(4.8) 


The critical wave number is determined by minimizing R (a) with 

C 

respect to a. Thus the critical wave number satisfies the equation 

T = 2 ( ^ (4.9) 


Thus introducing (4.9) in (4.8) determines the critical Rayleigh 
number . 


Re = -H- C (4.10) 

From equation (4.10) it is obvious that the minimum value of R 

c 

will be obtained for a = 0. This corresponds to R = 360. In the 

c o 

absence of rotation, we have 


R = 
c 


n 


8 


A rr2 ^ 2-2 

(a + a^ ) 


The minimum is once again obtained for a = 0 corresponding to 
which the R^=120. This is in complete agreement with the results 

1 7 

of Jakeman for stress free surfaces. Thus in the presence of 
rotation, the critical Rayleigh number is raised by a factor of 3 
if the fluid is confined within insulating boundaries 

.Corresponding to a^= 0 we have , 


R, ( 0 ) = ( 1 + T^) / 8 
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where R = 

i 


R 


n 


and T = -V 

* rr^ 

n 


R ( 0 ) = 

1C 


3 n‘ 


8 


Thus R^(0) will be a minimum if R^(0) < This 
correponds to 2 . For - 2 , we- therefore- have 
stationary convection in the form of long wavelength rolls- For 
T >2 stationary convection takes place in the form of finite wave 

i 

length rolls. This threshold that we determine is independent of 
the Prandtl number . 

Next we determine the critical Rayleigh number for the onset 
of instability as oscillatory convection. We substitute p = iw 
which represents the marginal state and obtain the characteristic 
equation of the matrix . 


-iw^-oo^(2<y + 


„2 ^ 2 
n + a 


) + i 




( 1 + t ) + 


2 O’ a 

2 ^ 2 
n + a 


- 8 


r O' 1 

n2 J 


2 2,., ^2. 
o-^ a ( 1 + t ) 

2 2 
n a 


2 

o-^ 8 r 


0 


(4.11) 


n 


Equating the real and imaginary parts to zero we obtain 


( 1 + t2 ) + 

n^r zT 


, 2 . 2 .2 
( ^ + a ) 


a 


a 


2 . 2 

n a 


7T 


4 

(4.12) 


os 


R 


os 


2 

a 

, 2 . 2 
( rt -I- a 


R representing the 

OS 


where r 


with 
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4 . 6 Ruppers Lortz instability 

The stability of the long wave length rolls is determined by 
investigating the Ruppers Lortz instability of such rolls. 
Accordingly, we superimpose perturbations in the form of a set of 
rolls with the axis along an arbitrary direction. Thus we 
represent the set of rolls along an arbitrary direction with 
modes of the form cos (k x + k_y) for the ,x and y dependence. Here 

k^ and k^ are the components of a along the x and y directions 
respectively. 

A minimal model representation for w, C and & consistent with 
the conditions on the boundaries for being stress free and 
thermally insulating is 

w = sin nz a(t) cos ax + a^(t) cos (k^x + k^y) 

C = cos nz f(t) cos ax + f^(t) cos(k^x + k^y) j- (4.14b) 

& = b(t) cos ax + c(t) cos riz + b^(t) cos (k^^x + k^y) 

+ cos Ttz { c (t) COS (k X + k y + a x)} 

JL X ^ 

+ COS nz { c^Ct) cos (k^x + k^y - a x)} (4.14c) 



Introducing these representations in equations (4 . la)-(4 . Ic) , 
we obtain the following dynamical system, with the time rescaled 
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Y = 


li 

X 

- X + 

Y + 

t G ) 

(4.15a) 

2 

8 r X 

4_ 

cos 0 

(X^Z^ - 

2 2 ^ ^ 
+ a"^) 


— . -I_ 

2 

, c sin 0 

( - 

fT- 7 

) - X Z 

(4.15b) 

2 



Z = 


Z 


2 2 
+ a^) 


+ X Y + 


(4.15c) 


G = - a i G+tX) 


(4.15d) 


X = O' ( - X + Y + t G ) 

1 111 


(4.15e) 


s Y + 

,2 2 . 1 

(n + a ) 


8 r X Q 

_j_ + ■ qg . | _.g_ (X z^- X z^) 

7T 


+ ■= ^ G ( Z^- 2^ ) 




(4-15f) 


z. = 


- Z - cos 0 (X Y + X Y ) + c sin 0 (Y G - Y G ) 

+ 1 1 1 . IX 

(4.15g) 


= - b_ Z^ + cos 0 (X^Y + X Y^ ) +c sin 0 (Y^G - Y G^ ) 

(4.15h) 

The rolls along the y axis are described by X,Y,Z,G non 
zero and the others zero. We consider the other set of rolls as the 
perturbation . A study of tine evolution of these perturbations 
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indicate that the condition for the onset of Khppers Lortz 
instability is given as 

(c t sin & - cos 0 ) ( - 2^) = 0 (4.16) 

c t sin 0 = cos 0 is the condition for Kuppers Lortz 

instability .The minimum value of t thus obtained is zero and it 
corresponds to a critical angle of 90°. for t > 0, the basic 

pattern is destabilized by those rolls which are at an angle for 
which ( fig. 7 ) c t :> cot 0 . Finally for very high values of 
t the original roll is unstable to a large range of rolls. 
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APPENDIX A 


In this appendix, we reduce the hydrodynamic equations to a 
scalar form. The equation of motion for the convecting fluid in the 
presence of rotation is given by equations (2.5a)-(2.5b) . 

Taking the curl on both sides of equation (2.5a) we get 


+ {V X (v.V)v } = V X 

Where w is the vorticity 


(g oi 0) + V 


CO 


+ 20 


a V 

3 Z 


(A.l) 


We consider the z component of this equation and arrive at the 

expression for C ( 2.16b). 

2 3 3 xa 

< |-t ) C = - T 


3 z 


Considering the curl of A.l once again and taking its z-component 
results in the expression for w which is given by equation (2.16a). 

Thus the equations has been transformed into a scalar form 
which are simpler to handle. 
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APPEMDIX B 

i ) . Solvability criterion 
We consider an inhomogeneous equation 

L I f > = I h > (Bl) 

which is to be solved under the constraint 

L I g > = 0 (B2) 

A vector < g^ | can be constructed with the property 

< gj I L = 0 (B3) 

Then from equation ( Bl ) we have the solvability condition 
<gj|L|f> = <gj^|h> = 0 (B4) 

Thus the condition is equivalent to stating that the 

inhomogeneous equations can be solved if the inhomogeneity is 
orthogonal to the solutions of the adjoint homogeneous problem. 

ii ) Newton Raphson method 

We employ the Newton Raphson method for determining the 

critical T for the onset of Ruppers Lorts instability. This 

is an iteration method for solving the equations f(x) =0. A series 
is generated which gives a sequence x^, x^ . -x^ recursively 

(B5) 

using the relation ^ i ^ 

where f'(x, ) denotes the derivative of f(x) at x = x^. An initial 
guess x^ however is required to get the iteration started. 
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APPENDIX C 

Integration of ordinary differential equations 

He employ the Runge- Kutta method of order four for the 
numerical integration of the dynamical system which in general can 
be represented as 

= f(x,y) (Cl) 

The formulae for the step by step integration which finally 
leads to a time series are 

y , = y + — ^ + 2 k„ + 2 ko + k. ) (C2) 

n + l 'n 6 1 Z O 4 

X ^ = X + h 

ri + 1 n 

where h is the time interval per step. 


k- = h f ( X , y ) (C4) 

JL n n n 

ko = h f (x + 0.5 h , y + 0.5 k., ) (C5) 

2 n n n n 1 

k„ = h f (x + 0.5 h , y„ + 0.5 k.,) (C6) 

3 r. n n n Z 

k. = h f (x + h , y + k„ ) (C7) 

4 n n 3 


4 

the truncation error is estimated to be of the order of h . 
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APPENDIX - D 


Kuppers Lortz instability with rigid boundaries 


In this appendix, we study the effect of rigid boundaries on 
the onset of Kuppers Lortz instability. The boundary conditions to 
be satisfied by the velocity and vorticity nodes are 

d w 

" I 

(Dl) 


C = 


d z 
0 


at z = ± 


While the assunption of the boundaries being perfect thermal 
conductors leads to 

0=0 atz=±| 

The rolls along the y-axis are represented as 


w 


a(t) cos a X ( 


2 


z 


1 )2 
4 ^ 


(D2) 


C = f(t) cos axz(z^-j) 


(D3) 


0 = b(t) cos ax( z^-^)+ c(t) z ( z^ - ^ ) 


(D4) 


We substitute these expansions in the hydrodynamic equations 
and obtain a set of Lorenz equations which are given by 


X = o-C-X + Y+ tG) 


P, Y = 


- Y + r X - 


X Z 
28 


(D5) 

(D6) 


p Z = 

^3 


X Y - b Z 


(D7) 
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Pi 


G = - O' 


{Xt+(42 + a2)(j^ 


where r = 


189 ^ ^ r 2 _ 

«8 ( 10 + a 2> * ^ 1 504 ) 


(D 8 ) 


(D9) 


(a ^ + 24 a ^ + 504) 
(12 + a 


P^ = 


(10 + a 


P 3 = ( 10 +a , 

^2 ^ 12 

( a + 24 a ^ +504 ) 

Performing the linear stability analysis 

state (X=Y=Z=G=0), we determine the 

number for the onset of convection as 


b = 42 (10 + a 


s.bout the conduct ion 
*^^itical Rayleigh 


r 


1 + 


( 42 + a 


Thus we find that the threshold Rayleigh n»».L. 

number is reduced by 

the introduction of rigid boundaries { r = ( 1 . ^ 2 . 

+ t ) for free 

boundaries }. Hinimizing with respect to ' 

results in the 

determination of the critical wave number. 


( 42 + a ^ )^ ( 2 a ® + 34 a ^ - 5040 ) 
( a ^ + 20 a ^ + 420 ) 


12 (DIO) 


In order to study the KUppers Lortz instability, of these 
rolls, we need to introduce a y dependence 

+n the velocity, 

vorticity and temperature modes. Assuming 

= a cos e , = a ^ 

we now superimpose a set of rolls whose ^ 

*ijcis IS along an 

arbitrary direction on the original set. Considering modes upto 
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second order we have 

w = < 2^ - I a(t) cos a X + a^(t) cos ( k^x + k^Y ) + 

z a (t) cos ( kx + ky + ax ) + za(t) cos ( k^x + k^y - ax) 

2 12^ -J 

(Dll) 


C - ( ^2 _ 1 ) j- ^ cos a X + z f,(t) cos ( k^x + k^y ) + 

cos ( k^x + k^y + a x ) + fgC^) ^ ''' a x ) j 

(D12) 


e = [ ( - 1 ) I b(t) cos a X + b/t) ooa ( k.x + ) } + 

z < - i ) [ o<t) + c,<t) cos < kjX + + a > + 

o^(t) cos ( k,x + k^y - a X) } ] (D13) 

Substitution of these expansions in the hydrodynamic 
equations results in 


{ 


44+2a^( 1+ cos & ) 
( 12 + a ^ ) 


} = {- 


t G + Z ( 1 + cos 9 ) 
2 1 


M+ 


^ \ 


{ 


44 + 2 a 


( 1 - cos ® ) 


(12 


2 X 

a ) 


} 


t G + ^ 2 ^ 

- T *3 
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P G 
^1 2 


[ t + g; { 


^ltX^+G^^ 10 + 2 a ( 1 + cos S ) 


}] 


(D14c) 


P^G^ = -<>'j^tXg + Gg-^10 + 2a^(l - cos ^ ^ ^ j 


(D14d) 


P 2 Y, = r X^- - 


28 


X (1-2 cos 0) 


4. 5 ^ ^ ^ ^ 5 ^ — 5^ 'V + f'Y + 

+ 168 ^ ^2 ) + c Y + 


168 

c Y sin 0 
84 


X 23 ( 1+2 COS 0 ) 
168 


{ 


G 


} 


( 1-cos 0) (1+COS 0) 

(D14e) 


( 1+2 cos 0) 


c sin 0 


P, Z, = - b_ z, + (X Y^+ X^ Y) + - (G y^-Y G^ ) 


2 2 


+ r Xg (10+a 


(D14f ) 


P, Z = - b, Z^ + ^ ;g9S -- g) (X Y + X^ Y) + ^ - j-— (GY - Y G. ) 

2 1+1 O 11 O 11 


+ r X^-i|f- (10+a 2) 


(D14g) 


X = cy(-X + Y+ tG) 


(D14h) 


{ 


p G = - c/^tX + G(42 + a‘') 

1 * 


} 


(D14i) 


^ XZ v_Y 7 (1-2 cos 0 ) _ 
P 2 * “ ^ ^ M * ■“‘1 168 


X. Z2 


(1-2 COS0 ) 
168 


G^ c sin 0 _ ^ Y c sin 0 


168 


^ ^2 ^ 84 


{ 


(1- cos 0) (1 


I 

+ cos 0) J 


(D14J) 
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P, Z = X Y + X Y - b Z 

^ jL 1 


(D14k) 


<5- ( - X + Y + t ) 
11 1 


(D141) 


Pi ®i = - P 


{ *1 


+ ( 42 + a G 




(D14in) 


where = 3960 +|2a^( 1+ cos © ) J + 176 a ^ (1+cos 9) 


N = [ 


a ^ + 24 a ^ + 504 




4 2 1 

a ^ + 24 a 504 


] 


.. = { 


42+2a^(l±cos© ) l|lO+a^ 


} { 10 . a 2 } 


The condition for Kdppers Lortz instability is obtained as 


X Z 

168 




1-2 cos 9 + 


1+2 cos 9 - 


c 

(42 


t sin 9 1 

1 a i 


Y c sin © 
84 


. . Q 'V X Z f 

c t sxn ® 1 ^ I 

(42 H. a 2) J 188 t 


i \ 

\ (1 - COS ©) (1 + cos ® ) J 


Using equations (D14a)- (D14ii) we obtain the critical Taylor 

4 

number to be of the order of 23 rt . In the limit of large Prandtl 
number, thus the boundary conditions have almost negligible effect 
on the onset of Ktippers Lortz instability. The boundary conditions 
probably play a more important role in the limit of small Prandtl 


number . 
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